

thetas<-read.csv(file="thetas.csv")[,1:3]


pd<-paste(wd, "/Figures/Figures_Output/",sep="")

pdf(file=paste(pd,"thetas.pdf",sep=""), 12, 9)


par(mar=c(4,4.25,2,1))

plot(thetas[,1],thetas[,2],ylim=c(0,1),type="n",axes=F,xlab=expression("%"~ Delta ~ "GDP"),ylab="Pr(Retain Power)",cex.lab=1.5,xlim=c(-.5,.5))

grid()

axis(1,seq(-1,1,by=.1), labels=seq(-100,100,by=10))

axis(2,seq(0,1,by=.2),labels=seq(0,100,by=20))

lines(thetas[,1],thetas[,2],col="dark red",lwd=2)

lines(thetas[,1],thetas[,3],col="dark blue",lwd=2)

text(-.25,.2,"Democracies",cex=1.5)
text(-.3,.92,"Autocracies",cex=1.5)

#mtext("Average (De)Stabilizing Effect of Growth",3,cex=2.5)



data<-read.csv("~/Dropbox/learning_democracy/data/data_final.csv")
data<-data[data$year > 1950 & data$year < 2001,]

rug(data$y)

dev.off()
